Order parameter model for unstable multilane traffic flow 
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We discuss a phenomenological approach to the description of unstable vehicle motion on multilane 
highways that explains in a simple way the observed sequence of the phase transitions "free flow 
<-^ synchronized mode «-> jam" as well as the hysteresis in these transitions. We introduce a new 
variable called order parameter that accounts for possible correlations in the vehicle motion at 
different lanes. So, it is principally due to the "many-body" effects in the car interaction in contrast 
to such variables as the mean car density and velocity being actually the zeroth and first moments of 
the "one-particle" distribution function. Therefore, we regard the order parameter as an additional 
independent state variable of traffic flow. We assume that these correlations are due to a small 
group of "fast" drivers and taking into account the general properties of the driver behavior 
formulate a governing equation for the order parameter. In this context we analyze the instability 
of homogeneous traffic flow manifesting itself in the mentioned above phase transitions and giving 
rise to the hysteresis in both of them. Besides, the jam is characterized by the vehicle flows at 
different lanes being independent of one another. We specify a certain simplified model in order 
to study the general features of the car cluster self-formation under the phase transition "free fiow 
<-> synchronized motion" . In particular, we show that the main local parameters of the developed 
cluster are determined by the state characteristics of vehicle motion only. 



INTRODUCTION. MACROSCOPIC MODELS 
FOR MULTILANE TRAFFIC DYNAMICS 



The study of traffic flow formed actually a novel branch 
of physics since the pioneering works by Lighthill and 
Whitham J|I| , Richards 0, and, then, by Prigogine and 
Herman [HfT It is singled out by the fact that in spite 
of motivated, i.e. a non-physical individual behavior of 
moving vehicles (they make up a so-called ensemble of 
"self-driven particles", see, e.g., &fl), traffic flow ex- 
hibits a wide class of critical and self-organizationphe- 
nomena met in physical systems (for a review see [0-^ ) . 
Besides, the methods of statistical physics turn out to 
be a useful basis for the theoretical description of traffic 
dynamics [ p^ . 

The existence of a new basic phase in vehicle flow on 
multilane highways called the synchronized motion was 
recently discovered by Kerner and Rehborn pd| ] , impact- 
ing significantly the physics of traffics as a whole. In 
particular, it turns out that the spontaneous formation 
of moving jams on highways proceeds mainly through a 
sequence of two transitions: "free flow — > synchronized 
motion — + stop-and-go pattern" |12| . All the three traffic 
modes are phase states, meaning their ability to persist 
individually for a long time. Besides, the two transitions 
exhibit hysteresis [O-MJ, i.e., for example, the transi- 
tion from the free flow to the synchronized mode occurs 
at a higher density and lower velocity than the inverse 
one. As follows from the experimental data ]ll| , p^ , [l4[ 
the phase transition "free flow <-!■ synchronized mode" is 
essentially a multilane effect. Recently Kerner 0^^ as- 



sumed it to be caused by "Z"-like form of the overtaking 
probability depending on the car density. 

The synchronized mode is characterized by substan- 
tial correlations in the car motion along different lanes 
because of the lane changing maneuvers. So, to de- 
scribe such phenomena a multilane traffic theory is re- 
quired. There have been proposed several macroscopic 
models dealing with multilane traffic flow and based on 
the gas- kinetic theory p5|-p0[ , a compressible fluid model 
[21 1 generalizing the approach by Kerner and Konhauser 
[22,p3|, and actually a model [p5,p5| dealing with the 
time dependent Ginzburg-Landau equation. 

All these models describe traffic flow in terms of the 
car density p, mean velocity v, and, may be, the velocity 
variance 9 or ascribe these quantities to vehicle flow at 
each lane a individually. In other words, the quantities 
{p, w, 9}a are regarded as a complete system of the traffic 
flow state variables and if they are fixed then all the vehi- 
cle flow characteristics should be determined. The given 
models relate the self-organization phenomena actually 
to the vehicle flow instability caused by the delay in the 
driver response to changes in the motion of the nearest 
cars ahead. In fact, let us briefly consider their simplifled 
version (cf. [^,^) which, nevertheless, catches the basic 
features taken into account: 
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Here the former term on the right-hand side of Eq. 



a so-called pressure term, reflects dispersion effects due 
to the finite velocity variance of the vehicles and the 
latter one describes the relaxation of the current velocity 
within the time r' to a certain equilibrium value Li{p, 0}. 
In particular, for 

V{p,vM = pe~ri^, (3) 

where r; is a "viscosity" coefficient and the velocity vari- 
ance 9 is treated as a constant, we are confronted with the 
Kerner-Konhauser model |£^,^ . In the present form the 
relaxation time r' characterizes the acceleration capabil- 
ity of the mean vehicle as well as the delay in the driver 
control over the headway (see, e.g., p8|"P0|). The value 
of the relaxation time is typically estimated as r' ~ 30 s 
because it is mainly determined by the mean time of vehi- 
cle acceleration which is a slower process than the vehicle 
deceleration or the driver reaction to changes in the head- 
way. The equilibrium velocity W{p} (here the fixed veloc- 
ity variance 6 is not directly shown in the argument list) 
is chosen by drivers keeping in mind the safety, the readi- 
ness for risk, and the legal traffic regulations. For homo- 
geneous traffic fiow the equilibrium velocity l^{p} = i9(p) 
is regarded as a certain phenomenological function meet- 
ing the conditions: 



dp 



< and pi?(jo) ^0 as p ^ po , 



(4) 



where po is the upper limit of the vehicle density on the 
road. Since the drivers anticipate risky maneuvers of 
distant cars also, the dependence U{p} is nonlocal. In 
particular, it is reasonable to assume that the driver be- 
havior is mainly governed by the car density p at a certain 
distant "interaction point" Xa = x+L* rather than at the 
current one x, which gives rise to a new term in Eq. (0) 
basing on the gas-kinetic theory |l^,|l^. Here, for the 
sake of simplicity following |g^ , we take this effect into 
account expanding p{x + L*) and, then, U{p} into the 
Taylor series and write: 



U{p} - d{p) - vo 



L^dp_ 
p dx 



(5) 



where vq is a certain characteristic velocity of the ve- 
hicles. Then linearizing the obtained system of equa- 
tions with respect to the small perturbations 5p^ Sv ex 
exp(7i -|- ikx) we obtain that the long-wave instability 
wiU occur if (cf. ||,||,|6i) 



r' {p^'f > VoL* + ' 



(6) 



in the long- wave limit the instability increment Re 7 de- 
pends on k as 

Re7 = fc2 [r'ip^'pf - {v„L* + t'9)] , (7) 

and the upper boundary fcmax of the instability region in 
the /c-space is given by the expression: 

1/2 
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(here -d' = d-d{p)/dp). As follows from (|^) the instabil- 
ity can occur when the delay time r' exceeds a certain 
critical value Tc and for r' < Tc the homogeneous traffic 
fiow is stable at least with respect to small perturbations. 
Moreover, the instability increment attains its maximum 
at fc ^ fcmax, so, special conditions are required for a wide 
vehicle cluster to form p2| , p3|j3jj34 | . In particular, in the 
formal limit r' — > from Eq. (0) we get 



^(p) 
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p dx 



(8) 



where D — vqL* plays the role of the diffusion coefficient 
of vehicle perturbations. The substitution of (||) into (||) 
yields the Burgers equation 



dp ^ d[pd{p)] _^d^p 

dt dx dx"^ 



(9) 



which describes the vehicle fiux stable, at least, with re- 
spect to small perturbations in the vehicle density p. 

However, the recent experimental data pi|-p4 33 
about traffic flow on German highways have demon- 
strated that the characteristics of the multilane vehicle 
motion are more complex (for a review see |R-§|). In 
particular, there are actually three types of synchronized 
mode, the totally homogeneous state, homogeneous-in- 
speed and totally heterogeneous ones flit]. Especially 
the homogeneous-in-speed state demonstrates the fact 
that, in contrast to the free flow, there is no direct rela- 
tionship between the density and flux of vehicles in the 
synchronized mode, because their variations are totally 
uncorrelated [[ll| . For example, an increase in the ve- 
hicle density can be accompanied by either an increase 
or decrease in the vehicle flux, with the car velocity be- 
ing practically constant. As a result, the synchronized 
mode corresponds to a two-dimensional region in the 
flow-density plane (jp-plane) rather than to a certain 
line j = i}(p)p pi| . Keeping in mind a hypothesis by 
Kerner [p|j9|j34[] about the metastability of each particu- 
lar state in this synchronized mode region it is natural to 
assume that there should be at least one additional state 
variable affecting the vehicle flux. The other important 
feature of the synchronized mode is the key role of some 
cars bunched together and traveling much faster than the 
typical ones, which enables us to regard them as a special 
car group [ pl| . Therefore, in the synchronized mode the 
function of car distribution in the velocity space should 
have two maxima and we will call such fast car groups 
platoons in speed. 

Anomalous properties of the synchronized mode have 
been substantiated also in | |33| using single-car-data. In 
particular, as the car density comes to the critical value 
Pc of the transition "free flow <-^ synchronized mode" the 
time-headway distribution exhibits a short time peak (at 
0.8 sec). This short time-headway corresponds to "...pla- 
toons of some vehicles traveling very fast-their drivers 
are taking the risk of driving "bumper-to-bumper" with 
a rather high speed. These platoons are the reason for 



occurrence of high-flow states in free traffic" |33|. The 
platoons are metastable and their destruction gives rise 
to the congested vehicle motion |35 . In the synchro- 



nized mode the weight of the short time-headways is less, 
however, almost every fourth driver falls below the 1-sec- 
threshold. In the vicinity of the transition "free flow 
<-> synchronized mode" the short time-headways have 
the greatest weight. In other words, at least near the 
given phase transition the traffic flow state is to be char- 
acterized by two different driver groups which separate 
from each other in the velocity space and, consequently, 
in multilane traffic flow there should be another relax- 
ation process distinct from one taken into account by the 
model (El), (0). In order to move faster than the sta- 
tistically averaged car a driver should permanently ma- 
noeuvring pass by the cars moving ahead. Meeting of 
several such "fast" drivers seems to cause the platoon 
formation. Obviously, to drive in such a manner requires 
additional efforts, so, each driver needs a certain time r 
to get the decision whether or not to take part in these 
manoeuvres. Exactly the time r characterizes the relax- 
ation processes in the platoon evolution. It should be 
noted that the overtaking manoeuvres are not caused by 
the control over the headway distance and, thus, the cor- 
responding transient processes may be much slower then 
the driver response to variations in the headway to pre- 
vent possible traffic accidents. 

The analysis of the obtained optimal-velocity func- 
tion V{Ax) demonstrates its dependence not only on the 
headway Aa; but also on the local car density. So, in con- 
gested flow the drivers supervises the vehicle arrangement 
or, at least, try to do this in a sufficiently large neighbor- 
hood covering several lanes. 

Another unexpected fact is that the synchronized 
mode is mainly distinctive not due to the car velocities 
at different lanes being equal. In the observed traffic 
flow various lanes did not exhibit a substantial difference 
in the car velocity even in the free flow. In agreement 
with the results obtained by Kerner 1 11 1 the synchronized 
mode is singled out by small correlations between fluctu- 
ations in the car flow, velocity and density. There is only 
a strong correlation between the velocities at different 
lanes taken at the same time, however, it decreases suffi- 
ciently fast as the time difference increases. By contrast, 
there are strong long-time correlations between the ffow 
and density for the free ffow as well as the stop-and-go 
mode. In these phases the vehicle flow directly depends 
on the density. 

Thereby, the free flow, the synchronized mode and the 
jammed motion seem to be qualitatively distinct from 
one another at the microscopic level. So, it is likely that 
to describe macroscopically traffic phase transitions the 
set of the state variables {p, v, 9}a should be completed 
with an additional parameter (or parameters) reflecting 
the internal correlations in the car dynamics. In other 
words, this parameter has to be due to the "many-body" 
effects in the car interaction in contrast to such external 
variables as the mean car density and velocity being ac- 
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FIG. 1. Transition region separating, e.g., the free-flow and 
synchronized mode. 



tually the zeroth and first moments of the "one-particle" 
distribution function. Thus, it can be regarded as an in- 
dependent state variable of traffic flow. The derivation 
of macroscopic traffic equations based on a Boltzmann 
like kinetic approach |Q has also shown that there is 
an additional internal degrees of freedom in the vehicle 
dynamics. 

In any case a theory of unstable traffic flow has to 
answer, in particular, to a question of why its two phases, 
e.g., the free flow and the synchronized mode, can coexist 
and, thus, what is the difference between them as well as 
why the separating transition region (Fig. nl) does not 
widen but keeps a certain thickness. Besides, it should 
specify the velocity u of this region depending on the 
traffic phase characteristics. There is a general expression 
relating the transition region velocity u to the density and 
mean velocity of cars in the free flow and a developed car 
cluster: p^, uj, and pd-, Vd, respectively, that follows 
from the vehicle conservation |l|, namely, the Lighthill- 
Whitham formula: 



PrjVcl ~ PfVf 
Pel - Pf 



(10) 



A specific model is to give additional relationships be- 
tween the quantities u, p/, Vf, and pd, Vd resulting from 
particular details of the car interaction. We note that 
a description similar to Eqs. (P, (||) dealing solely with 
the external parameters {/?, f } do not actually make a 
distinction between the free-flow and congested phases 
and their coexistence is due to the particular details of 
the car interaction. 

The transition "free flow <-> synchronized motion" is 
rather similar to aggregation phenomena in physical sys- 
tems such as undercooled liquid when in a metastable 
phase (undercooled liquid) the transition to a new or- 
dered (crystalline) phase goes through the formation of 
small clusters. Keeping in mind this analogy Mahnke et 
al. p7|-|39| have proposed a kinetic approach based on 
stochastic master equation describing the synchronized 
mode formation that deals with individual free cars and 
their clusters. The cluster evolution is governed by the 
attachment and evaporation of the individual cars and 
the synchronized mode is regarded as the motion of a 
large cluster. 

To describe such phenomena in physical systems it was 
developed an effective macroscopic approach, called the 
Landau phase transition theory M], that introduces a 
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FIG. 2. Schematic illustration of the car arrangement in 
the various phases of traffic flow and the multilane vehicle 
interaction caused by car overtaking. 



certain order parameter h characterizing the correlations, 
e.g., in the atom arrangement. In the present paper fol- 
lowing practically the spirit of the Landau theory we de- 
velop a phenomenological approach to the description of 
the traffic flow instability that ascribes to the vehicle flux 
an additional internal parameter which will be also called 
the order parameter h and allows for the effect of lane 
changing on the vehicle motion. In this way the free flow 
and the congested phases become in fact distinctive and 
solely the conditions of their coexistence as well the dy- 
namics of the transition layer are the subject of specific 
models. 



II. ORDER PARAMETER AND THE 
INDIVIDUAL DRIVER BEHAVIOR 

We describe the vehicle flow on a multilane highway in 
terms of its characteristics averaged over the road cross- 
section, namely, by the car density p, the mean velocity 
w, and the order parameter h. The latter is the measure 
of the correlations in the car motion or, what is equiv- 
alent, of the car arrangement regularity forming due to 
the lane change by the "fast" drivers. Let us discuss the 
physical meaning of the order parameter h in detail con- 
sidering the free flow, synchronized mode and jammed 
traffic individually (Fig. 0). 



A. Physical meaning of the order parameter h and 
its governing equation 

When vehicles move on a multilane highway without 
changing the lanes they interact practically with the 
nearest neighbors ahead only and, so, there should be 
no internal correlations in the vehicle flow at different 
lanes. Therefore, although under this condition the traf- 
fic flow can exhibit complex behaviour, for example, the 
"stop-and-go" waves can develop, it is actually of one- 
dimensional nature. In particular, the drivers that would 



prefer to move faster than then the statistically mean 
driver will bunch up forming the platoons headed by a 
relatively slower vehicle. When the cars begin to change 
lanes for overtaking slow vehicles the car ensembles at 
different lanes will affect one another. The case of this 
interaction is due to that a car during a lane change ma- 
noeuvre occupies, in a certain sense, two lanes simulta- 
neously, affecting the cars moving behind it at both the 
lanes. Figure Ep illustrates this interaction for cars 1 and 
2 through car 4 changing the lanes. The drivers of both 
cars 1 and 2 have to regard car 4 as the nearest neigh- 
bor and, so, their motion will be correlated during the 
given manoeuvre and after it during the relaxation time 
t'. In the same way car 1 is affected by car 3 because 
the motion of car 4 directly depends on the behavior of 
car 3. The more frequently lane changing is performed, 
the more correlated traffic flow on a multilane highway. 
Therefore, it is reasonable to introduce the order param- 
eter h being the mean density of such car triplets normal- 
ized to its maximum possible for the given highway and 
to regard it as a measure of the mulitlane correlations in 
the vehicle flow. 

On the other hand the order parameter h introduced 
in this way can be regarded as a measure of the vehi- 
cle arrangement regularity. Let us consider this ques- 
tion in detail for the free flow, synchronized mode, and 
jammed traffic individually. In the free flow the feasi- 
bility of overtaking makes the vehicle arrangement more 
regular because of platoon dissipation. So as the order 
parameter h grows the free traffic becomes more regular. 
Nevertheless, in this case the density of the car mulitlane 
triplets remains relatively low, /i <C 1, and the vehicle 
ensembles should exhibit weak correlations. Whence it 
follows also that the mean car velocity d is an increas- 
ing function of the order parameter h in the free flow. 
In the jammed motion (Fig. gc) leaving current lanes is 
hampered because of lack of room for the manoeuvres. 
So the car ensembles at different lanes can be mutually 
independent in spite of individual complex behavior. In 
the given case the order parameter must be small too, 
/i <C 1, but, in contrast, the car mean velocity should 
be a decreasing function of h. In fact, for highly dense 
traffic any lane change of a car requires practically that 
the neighbor drivers decelerating give place for this ma- 
noeuvres. 

Figure S illustrates the transition "free flow <-> syn- 
chronized mode" . As the car density grows in free flow, 
the "fast" drivers that at first overtake slow vehicles in- 
dividually begin to gather into platoons headed by more 
"slow" cars among them but, nevertheless, moving faster 
than the statistically mean vehicle (Fig. ga). The pla- 
toons are formed by drivers preferring to move as fast as 
possible keeping short headways without lane changing. 
Such a state of the traffic flow should be sufficiently in- 
homogeneous and the vehicle headway distribution has 
to contain a short headway spike as observed experimen- 
tally in p^]. Therefore, even at a sufficiently high car 
density the free flow should be characterized by weak 
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FIG. 3. Schematic illustration of the alteration in the vehi- 
cle arrangement near the phase transition "free flow ^-> syn- 
chronized mode". 



iTLultilane correlations and not too great values of the or- 
der parameter hf. The structure of these platoons is also 
inhomogeneous, they comprise cars whose drivers would 
prefer to move at different headways (for a fixed velocity) 
under comfortable conditions, i.e., when the cars moving 
behind a given car do not jam it or none of the vehicles 
moving on the neighboring lanes hinders its motion at the 
given velocity provided it changes the current lane. So, 
when the density of vehicles attains sufhciently high val- 
ues and their mean velocity decreases remarkably with 
respect to the velocity on the empty highway some of 
the "fast" drivers can decide that there is no reason to 
move so slowly at such short headways requiring strain. 
Then they can either overtake the car heading the cur- 
rent platoon by changing lanes individually or leave the 
platoon and take vacant places (Fig. ^). The former 
has to increase the multilane correlations and, in part, 
to decrease the mean vehicle velocity because the other 
drivers should give place for this manoeuvres in a suffi- 
ciently dense traffic flow. The latter also will decrease the 
mean vehicle velocity because these places were vacant 
from the standpoint of sufficiently "fast" drivers only but 
not from the point of view of the statistically mean ones 
preferring to keep longer headways in comparison with 
the platoon headways. Therefore, the statistically mean 
drivers have to decelerate, decreasing the mean vehicle 
velocity. The two manoeuvre types make the traffic flow 
more homogeneous dissipating the platoons and smooth- 
ing the headway distribution (Fig. Bp and the low frag- 
ment). Besides, the single- vehicle experimental data ||33|] 
show that the synchronized mode is singled out by long- 
distant correlations in the vehicle velocities, whereas the 
headway fluctuations are correlated only on small scales, 
which justifies the assumptions of the synchronized mode 



being a more homogeneous state than the free flow. We 
think that the given scenario describes the synchronized 
mode formation which must be characterized by a great 
value of the order parameter, hg > hf, and a lower veloc- 
ity in comparison with the free flow at the same vehicle 
density. 

In addition, whence it follows that, first, the left 
boundary of the headway distribution should be approx- 
imately the same for both the free flow and the synchro- 
nized mode near the phase transition, which corresponds 
the experimental data [ p3[ . Second, since in this case the 
transition from the free flow to the synchronized mode 
leads to the decrease in the mean velocity, the "fast" 
driver will see no reason to alter their behaviour and to 
move forming platoons again until the vehicle density 
decreases and the mean velocity grows enough. It is rea- 
sonable to relate this characteristics to the experimen- 
tally observed hysteresis in the transition "free flow <-> 
synchronized mode" |lj-|lj]. Third, for a car to be able 
to leave a given platoon the local vehicle arrangement at 
the neighboring lane should be of special form and when 
an event of the vehicle rearrangement occurs its follow- 
ing evolution depends also on the particular details of the 
neighboring car configuration exhibiting substantial fiuc- 
tuations. Therefore, the synchronized mode can comprise 
a great amount of local metastable states and correspond 
to a certain two-dimensional region on the flow-density 
plane (j'p-plane) rather than a line j — d{p)p, which 
matches the experimental data [ pT[ and the modern no- 
tion of the synchronized mode nature [0-H . This feature 
seems to be similar to that met in physical media with lo- 
cal order, for example, in glasses where phase transitions 
are characterized by wide range of controlling parame- 
ters (temperature, pressure, etc.) rather than their fixed 
values (see, e.g., pH). 

This uncertainty of the synchronized mode, at least 
qualitatively, may be regarded as an effect of the internal 
fiuctuations of the order parameter h and at the first step 
we will ignore them assuming the order parameter h to 
be determined in a unique fashion for fixed values of the 
vehicle density p and the mean velocity v. Thus for a 
uniform vehicle fiow we write: 



T— = -$(/i,p,w) 



(11) 



where r is the time required of drivers coming to the 
decision to begin or stop overtaking manoeuvres and the 
function ^(/i, p, v) possesses a single stationary point h = 
h{p,v) being stable and, thus. 



5$ ^ 



(12) 



The latter inequality is assumed to hold for all the val- 
ues of the order parameter for simplicity. We note that 
equation (|Il|) also allows for the delay in the driver re- 
sponse to changes on the road. However, in contrast with 
models similar to (|^) and (Q), here this effect is not the 



origin of the traffic fiow instability and, thus, its partic- 
ular description is not so crucial. Moreover, as discussed 
in the Introduction, the time r characterizes the delay in 
the driver decision concerning to the lane changing but 
not the control over the headway, enabling us to assume 
T> r'. 

The particular value h{v, p) of the order parameter re- 
sults from the compromise between the danger of acci- 
dent during changing lanes and the will to move as fast as 
possible. Obviously, the lower is the mean vehicle veloci- 
ty V for a fixed value of p, the weaker is the lane changing 
danger and the stronger is the will to move faster. Be- 
sides, the higher is the vehicle density p for a fixed value 
of V, the stronger is this danger (here the will has no 
effect at all). These statements enable us to regard the 
dependence h{v,p) as a decreasing function of both the 
variables v, p (Fig. ^ and taking into account inequal- 
ity (n2[) to write: 



^- >0, 7^ >0, 
av op 



(13) 



with the latter inequality stemming from the danger ef- 
fect only. 
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FIG. 4. Qualitative sketches of the order parameter h as 
a function of the vehicle mean velocity v and the density p 
specified by the behavior of individual drivers. 

Equation ( pd] ) describes actually the behavior of the 
drivers who prefer to move faster than the statistically 
mean vehicle and whose readiness for risk is greatest. Ex- 
actly this group of drivers govern the value of h. There is, 
however, another characteristics of the driver behavior, 
it is the mean velocity v — ^{h, p) chosen by the statisti- 
cally averaged driver taking into account also the danger 
resulting from the frequent lane changing by the "fast" 
drivers. This characteristics is actually the same as one 
discussed in the Introduction but depends also on the 
order parameter h. So, as a function of p it meets condi- 
tions (Q). Concerning the dependence of i?(/i, p) on h we 
can state that generally this function should be increas- 
ing for small values of the car density, p <^ po, because 
in the given case the lane changing practically makes no 
danger to traffic and all the drives can overtake vehicle 
moving at lower speed without risk. By contrast, when 
the vehicle density is sufficiently high, p '^ po, only the 
most "impatient" drivers permanently change the lanes 
for overtaking, making an additional danger to the most 
part of other drivers. Therefore, in this case the veloc- 
ity ^{h, p) has to decrease as the order parameter h in- 
creases. For certain intermediate values of the vehicle 




order parameter, h 1 

FIG. 5. A qualitative sketch of the mean vehicle velocity vs. 
the order parameter h for several fixed values of the vehicle 
density p. 



density, p k, pc, this dependence is to be weak. Fig. |^ 
shows the velocity ^(/i, p) as a function of h for differ- 
ent values of p, where, in addition, we assume the effect 
of order parameter h g (0, 1) near the boundary points 
weak and set 



s- 



at ft, = and h = 1 . 



(14) 



We will ignore the delay in the relaxation of the mean 
velocity to the equilibrium value v — i?(/i, p) because the 
corresponding delay time characterizes the driver control 
over the headway and should be short, as already dis- 
cussed above. Then the governing equation (O) for the 
order parameter h can be rewritten in the form: 



T-j^ = -(t>{h,p); 



dcf 



(/.(/i,p) ='<i>[ft,p,^(ft,p)]. (15) 



For the steady state uniform vehicle flow the solution of 
the equation (j){h, p) — specifies the dependence h{p) 
of the order parameter on the car density. Let us, now, 
study its properties and stability. 



B. Nonmonotony of the h{p) dependence and the 
traffic flow instability 

To study the local characteristics of the right-hand side 
of Eq. ^^) we analyze its partial derivatives 



d(j) _ 9$ a$ di} 

dh dh dv dh 

d(i) _ a$ a$ dd 

dp dp dv dp 



(16) 
(17) 



As mentioned above, the value of d^/dp is solely due to 
the danger during changing lanes, so this term can be 
ignored until the vehicle density p becomes sufficiently 
high. In other words, in a certain region p < ph ^ Po 
the derivative d(j)/dp ^ {d^ / dv){d'd / d p) < by virtue 
of (Q) and ([l3|). So, the local behavior of the func- 
tion h{p) (meeting the equality d0 = and, thus, 
dh/dp = —{d<f)/dp){d<f)/dh)~^) depends directly on the 
sign of the derivative dc/j/dh, it is increasing or decreasing 
for dcp/dh > or d(p/dh < 0, respectively. 




p. p. P, P;, 

vehicle density, p 
FIG. 6. The region of the traffic flow instabihty in the 
ftp-plane and the form of the curve h{p) displaying the de- 
pendence of the order parameter on the vehicle density. The 
plot is a qualitative sketch. 



For long- wave perturbations oc exp{ikx} of the car dis- 
tribution on a highway the density p can be treated as 
a constant at the lower order in k. Therefore, accord- 
ing to Eq. (nsf) the steady state traffic flow is unstable if 
d(j)/dh < 0. 

Due to dl^ ) and (|lj) the first term on the right-hand 
side of ( p^ is dominant in the vicinity of the lines h — 
and h = 1, thus, in this region the curve h{p) is increas- 
ing and the stationary state of the traffic flow is stable. 
For p < pc the value d-d/dh > (Fig. H), therefore, the 
whole region {0</i<l, 0<p< pc} corresponds to the 
stable car motion. However, for p > pc there can be 
a region of the order parameter h where the derivative 
d(f>/dh changes the sign and the vehicle motion becomes 
unstable. Indeed, the solution v = r]{h, p) of the equa- 
tion <i>(/i, p, w) — can be regarded as the mean vehicle 
velocity controlled by the "fast" drivers and is decreasing 
function of ft. because oidrj/dh = —{d^/dh)/{d^/dv)~^. 
So, once such "active" drivers become to change the lanes 
to move faster, they will do this as frequently as possi- 
ble especially if the mean velocity decreases, which corre- 
sponds to a considerable increase in h for a small decrease 
in V. So, it is quite natural to assume that the value of 
drj/dh for p > pc is sufHciently small and 



dh 



dv 



dd 
'dh 



drj 
'dh 



<0. 



(18) 



Under these conditions the instability region does exist, 
the curve h{p) can look like "S" (Fig. g) and its decreasing 
branch corresponds to the unstable vehicle flow. The 
lower increasing branch matches the free flow state of 
the car motion, whereas the upper one should be related 
to the synchronized mode because it is characterized by 
the order parameter coming to unity. 



C. Hysteresis and the fundamental diagram 

The obtained dependence h{p) actually describes the 
first order phase transition in the vehicle motion. Indeed, 
when increasing the car density exceeds the value pi the 
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FIG. 7. The mean vehicle velocity (o) and the vehicle flux 
(6) vs. the vehicle density for the limit values of the order pa- 
rameter /i = and /i = 1 as well as the resulting fundamental 
diagram (c). A qualitative sketch. 



free flow becomes absolutely unstable and the synchro- 
nized mode forms through a sharp jump of the order pa- 
rameter. If, however, after that the car density decreases 
the synchronized mode will persist until the car density 
attains the value p2 < pi- It is a typical hysteresis and 
the region {p2, Pi) corresponds to the metastable phases 
of traffic flow. 

Let us, now, discuss a possible form of the fundamen- 
tal diagram j = j{p) showing the vehicle flux j = p'd[p\ 
as a function of the car density p, where, by definition, 
■d[p\ = d[h{p),p]. It should be pointed out that here we 
confine our consideration to the region of not too large 
values of the car density, p < ph, where the transition 
"free flow <-> synchronized mode" takes place. The tran- 
sition "synchronized mode <-> jammed traffic" will be dis- 
cussed below. Fig. 0a displays the dependence d{h, p) of 
the mean vehicle velocity on the density p for the fixed 
limit values of the order parameter h = and 1. For a 
small values of p these curves practically coincide with 
each other. As the vehicle density p grows and until it 
comes close to the critical value Pc when the lane change 
danger becomes substantial, the velocity i?(l,p) practi- 
cally does not depend on p. So at the point pc at which 
the curves i9(l,/ci) and i?(0,p) meet each other, i9(l,p) is 
to exhibit sufficiently sharp decrease in comparison with 
the latter one. Therefore, on one hand, the function 
ji(p) = p^(l,p) has to be decreasing for p > pc- On 
the other hand, at the point pc for ft <C 1 the effect of 
the lane change danger is not extremely strong, it only 
makes the lane change ineffective, dd/dh ~ (Fig. ||). So 
it is reasonable to assume the function jo(p) = P'^i^^p) 
increasing near the point pc- Under the adopted assump- 
tions the relative arrangement of the curves jo{p), Ji(p) 
is demonstrated in Fig. M6, and Fig. Mc shows the funda- 
mental diagram of traffic flow resulting from Fig. m and 



Fig. 06. 

Concluding the present section we note that in the 
given description of the driver behavior governing the or- 
der parameter h the vehicle flux j{h,p) = p'd{h, p) is an 
external characteristics of traffic flow. So, the obtained 
form of the fundamental diagram does not follows di- 
rectly from the developed model, but can be interpreted 
sufficiently reasonable. It can be rigorously justified if the 
critical point pc corresponds to the maximum of the flux 
j{h* , p) for a certain fixed value h* of the order parame- 
ter. In other words, when the road capacity is exhausted 
and the following increase in the vehicle density leads to 
a decrease in the vehicle fiux the drivers divide into two 
groups, the majority prefer to move at their own lanes 
whereas the most "impatient" drivers change the lanes as 
frequently as possible, giving rise to the traffic instability. 
This problem, however, deserves an individual investiga- 
tion. 



III. PHASE COEXISTENCE. DIFFUSION 
LIMITED CLUSTER MOTION 

The previous section has considered uniform traffic 
fiow, so, analyzed actually the individual characteris- 
tics of the free flow and the synchronized mode. In the 
present section we study their coexistence, i.e., the condi- 
tions under which a car cluster of finite size forms. This 
problem, however, requires that the traffic flow model 
be defined concretely. Therefore, in what follows we 
will consider a certain simple model which illustrates the 
characteristic features of the car cluster self-organization 
without complex mathematical manipulations. 

As before, the model under consideration assumes the 
mean velocity relaxation to be immediate and modifies 
the governing equation dl3) in such a way as to ascribe 
the order parameter ft, to a local car group. In other 
words, we describe the vehicle fiow by the Lighthill- 
Whitham equation with dissipation (see, e.g., |42[ and 
also Introduction), replace the time derivative in Eq. ( p5| ) 
by the particle derivative, and take into account that 
the order parameter cannot exhibit substantial variations 
over scales / ~ 9^^^t < vqt {9 is the velocity variance, vq 
is the typical car velocity in the free flow). Namely, we 
write: 



dp ^ d[pd{h,p)] _^d^p 

dt dx dx^ 



dh 



d{h,p) 



dh 
dx 



C{h} 



(19) 
c^{Kp)+^{x,t). (20) 



Let us discuss the meaning of the particular terms of the 
given model. The Burgers equation (pj|), as already dis- 
cussed in Introduction, allows for the fact that drivers 
govern their motion taking into account not only the be- 
havior of the nearest cars, but the state of traffic flow 
inside the whole field of their front view of length. The ef- 
fective diffusivity D can be estimated as Z? ^ L*vq, where 



L* ^ Hs a front distance looked through by drivers as- 
sumed to be much greater that the scale Z, so 



Dt - IL* > f . 
The function (l){h, p) is of the form 



dcf 



cl){h,p) = h{l-h)[a{p)~h]^ 



(21) 



(22) 



where 



1 for p < Pc, 

a(p) ^ { (Pc + A - p)/A for Pc < p < Pc 







for p > pc + A . 



It describes such a driver behavior that h = and h = 1 
are the unique stable values of the order parameter for 
p < Pc and p > pc+A, respectively, whereas, for pc < p < 
Pc + A the points h = 0, h = 1 are both locally stable and 
there is an additional unstable stationary point, namely, 
h = a{p). The term 



dx^ ^dx 



(23) 



governs spatial variations in the field h{x,t) and takes 
into account that drivers mainly follow the behavior of 
cars in front of them and cars moving at the rear cannot 
essentially affect them. The mean car velocity depends 
on h and p as 

p^{h, p) ''^ p^o(l -h) + [pc-do - v{p - pc)]h . (24) 



The last term on the right-hand side of Eq. (^0|) charac- 
terizes the random fluctuations in the order parameter 
dynamics: 

(e(a;,i))=0, (25) 

{i{x,t)£,{x' ,t')) ^(j^lTS{x-x')S{t-t'), (26) 

where a is their dimensionless amplitude. Expressions 
( P^ ) and (24) gives the /i(yo)-dependence and the fun- 
damental diagram shown in Fig. simplifying the one 
presented in Fig 0. 
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FIG. 8. The dependence h{p) and the fundamental dia- 
gram of traffic flow described by the model ( |l9| ) , (120) . 

If we ignore the random fluctuations of the order pa- 
rameter h, i.e., set a = 0, then the model (p^), ( |20| ) will 
give us an artificially long delay (much greater than t) 
in the order parameter variations from, for example, the 



unstable point h — to the stable point h — I. Such 
a delay can lead to a meaningless great increase of the 
vehicle density in the free flow without phase transition 
to congestion. In order to avoid this artifact and to allow 
for the effect of real fluctuations in the driver behavior we 
also will assume the amplitude a to obey the condition 

iai: 



/ 



5/4 



<^<1 



(27) 



(cr ^ 1, because, otherwise, the traffic flow dynamics 
would be totally random). It should be noted that small 
random variations of the order parameter h near the 
points h = 0, h = I going into the regions h < and 
h > 1, respectively, do not come into conflict with its 
physical meaning as the measure of the car motion cor- 
relations. Indeed, the chosen values h = and h = 1 can 
describe a renormalization of real correlation coefficients 
h = hi>0,h2<l. 

According to Eq. (EQ), for the order parameter h the 
characteristic scale of its spatial variations is I, so, the 
layer 3^ separating the regions where /i sa and 1 is of 
thickness about I. Due to inequality ( pl| ) the car density 
on such scales can be treated as constant. Therefore, the 
transition region Cp between practically the uniform free 
flow and the congested phase is of thickness determined 
mainly by spatial variations of the vehicle density and on 
such scales the layer 3^ can be treated as an infinitely 
thin interface. In addition, the characteristic time scale of 
the layer Qh formation is about t, whereas it takes about 



the time r„ 



D/v'o 



t{L*/1) > T for the layer Cp 



to 



form. Thereby, when analyzing the motion of a wide car 
clusters we may regard the order parameter distribution 
h(x, t) as quasi-stationary for a fixed value of the car 
density p. 

Let us, now, consider two possible limits of the layer 
3/i motion under such conditions. 



A. Regular dynamics 

In the region pc < p < Pc + A until the value of a{p) 
comes close to the boundaries h = and h = 1 the effect 
of the random ffuctuations is ignorable. In this cas e by 
virtue of the adopted assumptions the solution of Eq. (pfl) 
that describes the layer S/^ moving at the speed u is of 
the form: 



h = 



1 



1 + tanh 



ut 



A 



(28) 



Here for the layer 9oi of the transition "free-ffow — > syn- 
chronized mode" and for the layer Sio of the opposite 
transition (Fig. ||) 



2^2 r- 

Aoi = / , Aio = ~2v2r]J 



(29) 
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FIG. 9. The distribution of the order parameter and the 
car density in the vicinity of the layers SJh of the transition 
between the free flow and the synchronized phase as well as 
the velocity of their motion vs. the local values pi of the car 
density. 



Uqi = ^0 - ^ y= [1 + 77^ - 2a{p^)] , 

^ V2r]vT 



"10 = ^0 



I 



2 V2t 
where we introduced the quantities 

A„ = z9(0,p,)-z9(l,p,), 



[2r]ya{pi) - {r]y - 1)] 



(30) 

(31) 



Vv 



^ \2V2i) 



,-| 1/2 



TAy 

2V21 



and Pi is the corresponding value of the car density inside 
the layers SJqi and ^ip- 

Expressions (|30|), (pl| ) describe the regular dynamics 
of the car cluster formation because the transition, for 
example, from the free ffow to the synchronized phase 
at a certain point x is induced by this transition at the 
nearest points. The dependence of the velocities uoi and 
Mio on the local car density pi is illustrated in Fig. H. The 
characteristic velocities attained in this type motion can 
be estimated as 



^n- 



max{?9oA/pc, 1/t} , 



so, under the adopted assumptions the regular dynamics 
does not allow for the sufficiently fast motion of the layers 
'ih upstream. 



B. Noise-induced dynamics 

As the car density p tends to the critical values pc or 
Pc + A the value of a{p) comes close to the boundaries 
a(pc) = 1 and a{pc + A) = 0, and the point h = I or 
h — becomes unstable, respectively. In this case the ef- 
fect of the random fluctuations ^(x, t) plays a substantial 



role. Namely, the phase transition, for example, from the 
free flow to the synchronized motion (for p « pc + A) is 
caused by the noise ^{x,t) and equiprobably takes place 
at every point of the region wherein p « pc + A rather 
than is localized near the current position of the layer 
5oi- Under these conditions the motion of the layers 
3/i can be qualitatively characterized by an extremely 
high velocity in both the directions, which is illustrated 
in Fig. H by dashed lines. 

We note that the noise-induced motion, in contrast to 
the regular dynamics, is to exhibit significant fluctuations 
in the displacement of the layer S/j as well as in its forms. 
This question is, however, a subject for an individual 
study. 



C. Diffusion limited motion of vehicle clusters 

Let us, now, analyze the motion of a sufficiently large 
cluster that can form on a highway when the initial car 
density or, what is the same, the average car density 
p belongs to the metastable region, p S {pc, Pc + A). 
The term "sufficiently large" means that the cluster di- 
mension L is assumed to be much greater than the front 
distance L* looked through by drivers, so, they cannot 
look round the congestion as a whole. Exactly, in this 
case a quasi-local description of traffic flow similar to the 
differential equations ([l9|), ( pOJ ) is justified. 

Converting to the frame y — x — ut moving at the clus- 
ter velocity m, solving Eq. (^ individually for the free 
flow and the synchronized phase, and treating the lay- 
ers 3h as inflnitely thin interfaces we get the following 
conclusion. Within the frameworks of the given model 
the car cluster moves upstream sufficiently fast, so, the 
motion of the layers 3oi f^nd 5io is governed by the noise 
^(x,t). In this case the values of the car density at the 
layers Qqi and Qio have to be pj « pc + A and pf ~ pc, 
respectively. Thereby, the cluster velocity u is mainly de- 
termined by the car redistribution governed by the diffu- 
sion type processes. The latter feature is the reason why 
we refer to the cluster dynamics under such conditions 
as to the diffusion limited motion. The transition region 
£01 between practically the uniform free flow state and 
the cluster contains the exponential increase of the vehi- 
cle density inside the free flow phase from the value p/ 
far from the "interface" Sqi up to pj « pc + A at 3oi j 

P = Pf + {Pj ~ Pf) exp{g/y} , 

where qf — {do + \u\)/D ^ 1/L* and the frame {y} is at- 
tached to the "interface" 3oi- The transition region £10 
from the synchronized phase to the uniform free flow is 
to be localized inside the car cluster. So, it is character- 
ized by the decrease in the vehicle density Sp ex expl^j-y}, 
where qj — {\u\ — v)/ D, and the vehicle free flow leaving 
the cluster is uniform at all its points (Fig. pl]|a). 

The cluster velocity is directly determined by the mo- 
tion of the interface Sqi- Therefore, assuming also the 



cluster dimension L large in comparison with L* , from 
Eq. (|9|) we get the expression of the same form as the 
Lighthill-Whitham formula ( p!o[ ) relating the cluster ve- 
locity u and the vehicle flux characteristics on the both 
sides of the layer £01 • Whence it follows that at the flrst 
approximation: 

uK-v, (32) 

the value qj = 0, and the vehicle cluster is of the form 
shown in Fig. |lO| a under the name "mesocluster" . Assum- 
ing the total number of cars on the highway of length Lj-d 
flxed we get the expression for the mesocluster dimension 
L: 



L = 2L 



P' 



rd- 



(33) 



However, this result is justified only for sufficiently 
small values of (p — Pc)/^ ^ 1 , when the cluster di- 
mension is not too large, Lqj <C 1 (nevertheless, L ^ 
L*). Exactly for this reason we refer to such clusters as 
mesoscopic ones. In order to study the opposite limit, 
Lqj 3> 1, we have to take into account that the value pf 
is not rigorously equal to pc but practically is the root 
pZ > Pc of the equation uio(pf ) = —v. In this case the 
Lighthill-Whitham formula (Uw gives the expression: 



.Pf 



(^0 



leading to the following estimates of the thickness 1/qj 
of the transition region £10: 



l/'Zj 



DA 



L* 



(1?0 + I^)(p}-Pc) ip}-Pc) 



The form of such a wide cluster is shown in Fig. |lO|a, its 
dimension is: 



L — Li-d 



P- Pc 



(34) 



and the region of the mean car density corresponding to 
this limit is specified by the inequality: 



p* -Pc L* A 



Lid {p*f - Pc) 



(35) 



The resulting dependence of the cluster dimension on 
the mean car density p is illustrated in Fig. plip . 




P _ p +A 
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FIG. 10. The possible forms of tlie car clusters and their 
dimension vs. the mean car density. 
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IV. PHASE TRANSITION "SYNCHRONIZED 
MODE ^ JAM". BRIEF DISCUSSION 



In Sec. O we have considered the phase transition be- 
tween the free flow and the synchronized mode. However, 
according to the experimental data |1^ there is an addi- 
tional phase transition in traffic flow regarded as the tran- 
sition between the synchronized motion and the jammed 
"stop-and-go" traffic. This transition occurs at extremely 
high vehicle densities p coming close to the limit value po • 

The present section briefly demonstrates that the de- 
veloped model for the driver behavior also predicts a sim- 
ilar phase transition at high car densities. To avoid possi- 
ble misunderstandings we, beforehand, point out that the 
model in its present form cannot describe details of the 
transition "synchronized mode <-*■ jam" because we have 
not taken into account the delay in the driver response to 
variations in headway. The latter is responsible for the 
formation of the "stop-and-go" pattern, so, to describe 
the jammed traffic on multilane highways we, at least, 
should combine a governing equation for the order pa- 
rameter h and a continuity equation similar to (EQ), (119) 
with an equation for the car velocity relaxation similar 
to (0). This question, however, is worthy of individual 
study. Besides, the approximations used in Sec. |lll| to 



Besides, the approximations used in Sec. 
characterize the synchronized mode at the car densities 
near p^ do not hold here. 

In Sec. O we have studied the dependence of the order 
parameter h on the car density ignoring the first term on 
the right-hand side of Exp. (O) caused by the danger- 
ous of lane changing. This assumption is justified when 
the car density is not to high. In extremely dense traf- 
fic flow, when the car density exceeds a certain value, 
P > Ph ^ POj changing lanes becomes sufficiently dan- 
gerous and the function ^{h,v,p) describing the driver 
behavior is to depend strongly on the vehicle density in 
this region. In addition, the vehicle motion becomes suf- 
ficiently slow. Under such conditions the former term on 
the right-hand side of expression (17) should be dominant 
and, thus, d<j)/dp > 0. Therefore, the stable vehicle mo- 
tion corresponding to d4>/dh > matches the decreasing 
dependence of the order parameter h(p) on the vehicle 
density p for p > ph- So, as the vehicle density p in- 
creases the curve h{p) can again go into the instability 
region (in the /ip-plane), which has to give rise to a jump 
from the synchronized mode with greater values of the 
order parameter to a new traffic state with its less val- 
ues (Fig. |ri|). Obviously, this transition between the two 
congested phases also exhibit the same hysteresis as one 
describe in Sec. O. 

We identify the latter traffic state with the jammed 
vehicle motion. Indeed, in extremely dense traffic lane 
change is practically depressed, making the car ensem- 
bles at different lanes independent of one another. So, 
in this case vehicle flow has to exhibit weak multilane 
correlations and we should ascribe to it small values of 
the order parameter h. It should be noted that the ex- 




vehicle density, p 
FIG. 11. The instability region and the /i(p)-dependence 
describing the transition from the synchronized (congested) 
phase to the heavy congested phase (a jam) in the region of 
high car density. 



perimental single-vehicle data |p3[ demonstrates strong 
correlations of variations in the traffic flux and the car 
density for both the free flow and the "stop-and-go" mo- 
tion. By contrast, the synchronized mode is character- 
ized by small values of the cross-covariance between flow, 
speed, and density. In other words, for the free flow and 
the "stop-and-go" motion the traffic flux j = dp should 
depends directly on the car density p, as it must in the 
present model if we set h — Q. 

Finalizing the present section we point out that the 
given model treats the jammed phase as a "faster" vehicle 
motion then the synchronized mode at the same values of 
the order parameter. There is no contradiction with the 
usual view on the synchronized mode as a high flux trafflc 
state. The latter corresponds to the traffic fiow at the 
vehicle densities near the phase transition "free flow <-> 
synchronized mode" rather than close to the limit value 
Po- Besides, ordinary driver's experience prompts that 
a highly dense traffic flow can be blocked at all if one 
of the cars begin to change the lanes. Nevertheless, in 
order to describe, at least qualitatively, the real features 
of the phase transition "synchronized mode <-> stop-and- 
go waves" a more sophisticated model is required. The 
present description only relates it to the instability of the 
order parameter at high values of the vehicle density. 

Besides, the present analysis demonstrates also the 
nonmonotonic behavior of the order parameter as the car 
density increases even if we ignore the hysteresis regions 
and focus our attention the stabel vehicle flow regions 
only. It should be noted that a similar nonmonotonic de- 
pendence of the lane change frequency on the car density 
as well as the platoon formation has been found in the 
cellular automaton model for two-lane traffic fp4| . 



V. CLOSING REMARKS 

Concluding the paper we recall the key points of the 
developed model. 

We have proposed an original macroscopic approach 
to the description of multilane traffic flow based on an 
extended collection of the traffic flow state variables. 
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Namely, in addition to such characteristics as the car 
density p and mean velocity v being actually the zeroth 
and first moments of the "one particle" distribution func- 
tion we introduce a new variable h called the "order pa- 
rameter" . It stands for the internal correlations in the 
car motion along different lanes that are due to the lane 
changing manoeuvres. The order parameter, in fact, al- 
lows for the essentially "many-body" effects in the car 
interaction so it is treated as an independent state vari- 
able. 

Taking into account the general properties of the driver 
behavior we have stated a governing equation for the 
order parameter. Based on current experimental data 
[jll|- p^ , p3[ we have assume the correlations in the car mo- 
tion on multilane highways to be due to a small group 
of "fast" drivers, i.e. the drivers who move substantially 
faster than the statistically mean vehicle continuously 
overtaking other cars. These "fast" cars, on one hand, 
increase individually the total rate of vehicle flow but, on 
the other hand, make the accident danger greater and, 
thus, cause the statistically mean driver to decrease the 
velocity. The competition of the two effects depends on 
the car density and the mean velocity and, as shown, 
can give rise to the traffic flow instability. It turns out 
that the resulting dependence of the order parameter on 
the car density describes in the same way the experimen- 
tally observed sequence of phase transitions "free flow ^^ 
synchronized motion <~> jam" typical for traffic flow on 
highways |jl2]. Besides, we have shown that both these 
transitions should be of the first order type and exhibit 
hysteresis, matching the experimental data |13-n4|. The 
synchronized mode is characterized by a large value of 
the order parameter, whereas the free flow and the jam 
match its small values. The latter feature enables us to 
treat the jam as a phase comprising the vehicle flows at 
different lane with weak mutual interaction because of 
the lane changing being depressed. 

In order to illustrate the characteristic features of the 
car clusters that self-organizing under these conditions 
we have considered a simple model dealing only with the 
evolution of the car density and the order parameter. In 
particular, it is shown that in the steady state the car 
density inside the cluster and the free flow being in equi- 
librium with the cluster, as well as the velocity at which 
the cluster moves upstream are fixed and determined by 
the basic properties of the traffic flow. On the contrary, 
the size of the car cluster depends on the initial condi- 
tions. 

Finally, we would like to underline that the developed 
model takes into account only one effect causing the traf- 
fic flow instability. The other, the delay in the driver 
control over the headway, seems to be responsible for 
the "stop-and-go" waves in the jammed phase (for a re- 
view of the continuum description of this phenomena see, 
e.g., [p4U23]). So, combining the two approaches into one 
model it enables detail description of a wide class of phe- 
nomena occurring in the transitions from free flow to the 
heavy congested phase on highways. In this way the or- 



der parameter model could describe also the formation 
of a local jam on a highway whose boundaries comprise 
both of the phase transitions. In the present form it fails 
to do this because the free flow and the jammed traffic 
are characterized by small values of the order parameter. 
Concerning a possible derivation of the order param- 
eter model from the gas-kinetic theory we note that the 
appearance of the "fast" driver platoons demonstrates 
a substantial deviation of the car distribution function 
from the monotonic quasi-equilibrium form. So, to con- 
struct an adequate system of equations dealing with the 
moments of the distribution function a more sophisti- 
cated approximation is required. 
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